% 采样点和权重, Gauss 4
pos = [
    1/3,1/3
    0.73333333,0.13333333
    0.13333333,0.73333333
    0.13333333,0.13333333];
coe = [
-0.56250000,...
0.52083333, 0.52083333, 0.52083333];